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Abstract 

A simple model of estimating the radial distribution function of an imperfect Bose 
gas in the ground state is presented. The model is based on integro-differential equations 
derived by considering the space boson distribution in an external field. With the approach 
proposed, the particular case of dilute Bose gas is investigated within the hard sphere 
approximation and beyond. 

PACS numbers: 05.30.Jp, 67.90. +z 



1 Introduction 

This letter is devoted to investigation of spatial particle correlations in many-body systems in 
the spirit of the approach of Ref.l which consists in the following. Let us consider a uniform 
many-body system with the interaction potential $(r) at the particle density n = N/V. As 
it has been shown for the classical case Jl|, the radial distribution function of this system g(r) 
can be calculated with 



where n str (r) denotes the density of the mentioned particles at the point f (r =| f |) of the 
nonuniform structure appearing under the action of the external field $(r) on the system 
considered. Here it is implied that 

$(r) — > , n str {r) — > n (r — > oo) . 

Relation ([[]) seems quite natural because the quantity g(\ fx — r 2 \)/V is the probability of 
finding the second particle of the system at the point r 2 provided that the first one is at the 
point fi. Therefore, we are able to calculate g{r) by considering the space distribution of the 
particles with the numbers 2, N in the external field generated by the first particle fixed at 
the origin of the coordinates. 

The aim of the present article is to clarify to what extent relation (|I]) is suitable for the 
case of an imperfect Bose gas at zero temperature. The paper is organized as follows. To 
simplify understanding of the subject, the results of investigating relation([TJ) in the classical 
case are briefly discussed in the second section. The third paragraph contains the main items of 
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deriving integro-differential equations for g(r) of a Bose gas in the ground state on the basis of 
([[]). To check validity of the equations, in the fourth section the Bose gas of the hard spheres is 
considered in the weak coupling approximation. At last, the fifth paragraph is the consideration 
of a dilute Bose gas beyond the hard sphere simplification. 



2 Classical Case 



To demonstrate reasonable character of relation (|I|) and simplify understanding of further 
arguments, let us consider a many-body system of classical particles interacting with each other 
by means of the potential $(r) and being in the external field $(r). In this case we are able to 
calculate the structure function n str {r) with the well-known Thomas-Fermi approximation ||. 
According to this method, n str {r) can be found with the following condition: 



H (n str (r)) + $(r) + U{r) = const (Vr) 



(2) 



where /i(n str (r)) is the chemical potential of the ideal gas of the given particles at the density 
^str(^) ; U(r) denotes the energy of the interaction of the particle at the point r with the 
surrounding particles. To complete the calculational procedure, we use the reasonable integral 
relation 

U(r)= [ <$>{\r-y\)n str {y)dy. (3) 



Now, we have got everything to obtain an integral equation for g(r) in the classical case. At 
some point ff ar being far enough from the origin of the external field, expression (|2]) is written 

as 

/j,(n) + n $(| f far - y \)dy = const 
J v 

provided that we limit ourselves to the situation 



(4) 



lim $(r) 

r— >oo 







Using formulae 



and the relation 



we arrive at 



the classical approximation 

fi(n) - fi(n str (r)) = -9 In (n str (r)/n) 



lim 



Jv(H\ r f ar - y\)-$(\r-y \))dy 
f v m f-y\)dy 



- 9 \ng(r) = $(r) + n / (g(y) - 1)$(| f-y \)dy 



v 



(5) 



(6) 



Note that relation (||) is valid for lots of the known potentials, in particular, for the integrable 
and Coulomb — like potentials. 

We have derived integral equation (||) neglecting spatial particle correlations when calculat- 
ing U(r) QTJ] . Therefore, it can only be used for the integrable and Coulomb-like potentials but 
not for the strongly singular potentials with the short range behaviour of the Lennard- Jones 
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type. To include the correlations, we should be based on a more accurate approximation of 
U(r). The simplest one is the following ]1||: 

U{r) = J g{\ f-y\) $(| f-y |) n str (y)dy . (7) 

With (D, (D and (0) we find 

- 9 lnp(r) = $(r) + n / - 1) p(| r - y |) $(| f - y |) . (8) 

In spite of the simplest version of approximating U(r) with the space particle correlations taken 
into account, integral equation (|8]) is very similar to the well-known Bogoliubov equation for 
g{r) (see discussion in Ref.l). As it is seen, relations (^|) and (||) can be used for the integrable, 
Coulomb-like as well as Lennard- Jones potentials because g(r)$(r) — > (r — > 0) for the 
strongly singular $(r). Thus, in this paragraph we have been convinced of the correctness of 
relation (p]) for the classical many-body systems. To elucidate the problem of using (p]) in the 
quantum case, we further investigate many-boson system in the ground state. 



3 Bose Gas in the Ground State 

One could expect that in the quantum case being essentially more complicated than the classical 
one, such a simple approach based on expression (|]) is not able to provide us with satisfactory 
results. However, in the situation with no exchange effects, for the Bose gas in the ground 
state, relation ([]]) yields a quite reasonable estimate of g(r) with only one obvious correction. 
To obtain more accurate data on g(r) for a uniform many-body system of bosons with the 
mass m, we need to investigate a nonuniform Bose gas made of particles with the mass m/2 . 
It should be noted, that as the quantity m does not appear in the integral equations of the 
previous paragraph, the mentioned correction does not contradict the classical case either. 

Now, let us consider the system of iV — 1 zero-spin bosons with the mass m/2 at zero 
temperature and in the external field $(r). As before, we assume the interparticle potential to 
be $(r). The simplest way to find n str (r) in this case is to adopt the Hartree approximation 
for the wave function of the system 

JV-l 

^ = n (9) 

i=i 

where ?/>(rj) denotes the normalized wave function of the i— th boson in the ground state. As 
ip(r) can be chosen to be a real function ||, we have 

n str {r) = (iV-l)^ 2 (r), (10) 

and, therefore, ip(f) = ip(r). Thus, to calculate n str (r), we need evaluating ip{r) which obeys 
the following Schrodinger equation: 

h 2 

Aip(r) + $(r)ip(r) + 

777* 

+(N-2)1>(r) f mr-y\)ij 2 (y)dy = E ^(r). (11) 

Jy 
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Here Eq denotes the energy of a boson in the ground state. Not to solve the eigenvalue problem, 
let us take into account that at the point r/ ar being far enough from the origin of the coordinates 

^{r far ) = ^= (*(r /or ) = 0). (12) 

Substituting ( |T2| ) into ( [XT] ) one can obtain 

iV-2 r 

~y- J v mr far -y\)dy = E . (13) 
Further, with the use of relations @, (|TT|) and ([H|) we arrive at 



h 2 



+ (^_^( r ) jf $(| r — y |) (V^(</) - l) dy. (14) 



Taking into account that in the thermodynamic limit 

N-2 N 



V V 



n 



and using ([!]), (|T0D and (0), we derive the following integro-differential equation for g(r) of a 
cold Bose gas: 

-/\g l l 2 {r) = $(r) + 



+n I $(| f-y (15) 



By analogy with the classical case this equation can be generalized to the situation of strongly 
singular potentials 

-Ag^ 2 (r) = $(r) + 



+n I g(\f-y |)$(| r - y |) - 1) dy. (16) 



v 



Remark that equations ([15|) and ( fig) have been derived for the spinless bosons but they are 
also valid for bosons with nonzero spin. 
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4 Bose Gas with the Hard Sphere Interaction 

To verify the reasonable character of equations (|T5| ) and ([16]), let us consider the well-known 
example of the hard-sphere Bose gas in the weak coupling approximation characterized by the 
relation 

g(r) = l — e(r), s(r) <C 1. 

Taking account of this relation and limiting oneself to the first order in e(r), one can rewrite 
equation fll~5|) (or (|l6l)) as 

h 2 f 
- —Ae(r) = $(r) -n $(| r - y \)e(y)dy. (17) 
2m Jv 

The new equation is solved with the Fourier transformation and gives the following result: 

I f $(o) exp(tqr 



(2tt) 3 J h 2 q 2 /2m + n$(q) 



dq 



where $(g) denotes the Fourier transform of $(r) . Inserting $(g) = 4:irh 2 a/m, where a denotes 
the sphere radius, we find 
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9{r) 



fest(x), x = 4r v /: 



nan 



(19) 



here f es t(x) = 0.5 exp(— xjy 2). How well this estimate is coordinated with the explicit expres- 
sion (see Ref.4) 

g (r) = l-^f expl (x), (20) 

is shown in the following table: 



X 


fest(%) 




X 


fest(x) 


fexpl 


X 


fest(%) 


fexpl (*^) 


0.0 


0.500 


0.500 


2.0 


0.122 


0.098 


5.0 


0.015 


0.012 


0.5 


0.351 


0.328 


2.5 


0.085 


0.067 


10 


0.0004 


0.001 


1.0 


0.246 


0.217 


3.0 


0.060 


0.046 


20 


3.6- 10" 6 


0.0002 


1.5 


0.173 


0.145 


4.0 


0.030 


0.023 









As it is seen, esti- 



mate (^) is in reasonable agreement with the explicit result (|20|) except the region of large 
distances 

10 



r > 



V2 



where the correlation length 



in an 



is the same for both the (19) and (|20|) expressions. Note that in the case of (19) the correlation 
length can be found from exp(— x/\/2) = exp(— 1) . Estimation (|T9"D does not give the known 
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1/r 4 —decay of g(r) for r>r c |||. Therefore, we need to be careful investigating the spectrum 
of the low- lying excitations in a Bose gas with g(r) estimated from equations ( Jl5| ) and (Jl6|). 
However, the model proposed can yield satisfactory evaluations for the main thermodynamic 
characteristics which are essentially determined by the values of g(r) at r < r c and r ~ r c . For 
example, the ground state energy calculated with (O) is given by 



2nh 2 a N 2 / a 3 N . 

E = — T 1 + 4 W— I- (21) 



The explicit result (BO) leads B to the expression 



2nh 2 a N 2 ( a 3 N\ , s 

e = — T + J (22) 

first found by Lee and Yang |7f . Besides, we can improve our estimates treating the boson mass 



in equations (O) and (|16D as a free parameter of the model. In particular, to derive (p2|) we 



should replace m by 1.13m in flT5j ) or 



5 Natural Dilute Bose Gas 



A helpful feature of equations (15) and (|16|) consists in the possibility of operating with them 
beyond the weak coupling and hard spheres. In particular, in investigating a natural dilute 
Bose gas the expansion in powers of n for g(r) may be of interest. To find it let us assume that 

gV\r) =u (r) +n Ul {r) + ... (23) 

and substitute (p3|) into, say, equation (]T6|) . We arrive at the following differential equations: 



and 



fly 2 . , , 

— Auur) 
m 



— Au (r) 
m 



$(r) u\{r) 



$(r) u (r) = 0, 



(24) 



Uo(r) / u 2 {\ r -2/|) $(| f - y |) (^(y) - l)dy. 



(25) 



Relation (24) is exactly the equation which is often used to find the zero density limit for the 
function of the boson two-body distribution 



Equations fl24| ) and ( p5[ ) testify to that g(r) 
of a natural Bose gas in the ground state is an analytical function of n contrary to the case of 
the weakly interacting Bose gas of the hard spheres. 
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6 Summary 



The integro-differential equations for the radial distribution function of a cold Bose gas with 
the interaction potential $(r), have been found via investigating the spatial boson distribution 
in the external field $(r) . Checking the validity of these equations, we have considered the 
weakly interacting Bose gas of the hard spheres and derived the reasonable picture of the two- 
body correlations and the satisfactory estimate of the mean energy. In the case of a natural 
Bose gas ( beyond the hard spheres and weak coupling ), the equations have provided us with 
the expansion in powers of n for g(r) that can be useful at investigating a dilute many-boson 
system. 

Acknowledgement: The author thanks A. Yu. Cherny for helpful and stimulating discussions. 
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